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ABSTRACT 

Milhmeter wavelength observations are presented of NGC 1333 IRAS 4, 
a group of highly-embedded young stellar objects in Perseus, that reveal 
motions of infall, outflow, rotation, and turbulence in the dense gas around its 
two brightest continuum objects, 4A and 4B. These data have finest angular 
resolution of ~2" (0.0034 pc) and finest velocity resolution of 0.13 km s^^. 

Infall motions are seen from inverse P-Cygni profiles observed in H2CO 
3i2-2ii toward both objects, but also in CS (3-2) and N2H+ 1-0 toward 4A, 
providing the least ambiguous evidence for such motions toward low-mass 
protostellar objects. Outflow motions are probed by bright line wings of H2CO 
3i2-2ii and CS (3-2) observed at positions offset from 4A and 4B, likely tracing 
dense cavity walls. Rotational motions of dense gas are traced by a systematic 
variation of the N2H"'" line velocities, and such variations are found around 4A 
but not around 4B. Turbulent motions appear reduced with scale, given N2H+ 
line widths around both 4A and 4B that are narrower by factors of 2 or 3 than 
those seen from single-dish observations. Minimum observed line widths of ~0.2 
km s^^ provide a new low, upper bound to the velocity dispersion of the parent 
core to IRAS 4, and demonstrate that turbulence within regions of clustered 
star formation can be reduced significantly. 

A third continuum object in the region, 4B', shows no detectable line 
emission in any of the observed molecular species. 



^currently at the Radio Astronomy Laboratory, 601 Campbell Hall, University of California, Berkeley, 
Berkeley, CA, 94705-3411 U.S.A. 

^currently at the Institute of Astronomy and Astrophysics, Academia Sinica, P.O. Box 1-87, Nankang, 
Taipei 11529, Taiwan, R.O.C. 

■^currently at the Dcpartamento de Astronomia, Universidad de Chile, Casilla 36-D, Santiago, Chile 



- 2 - 



Subject headings: ISM: individual (NGC 1333 IRAS 4), ISM: kinematics and 
dynamics, stars: formation 



-3- 



1. Introduction 

Star formation involves both coherent and incoherent gas motions within molecular 
clouds. Coherent motions are those with an implicitly common directionality, and include 
the infall, outflow, and rotational motions of molecular gas associated with extremely 
young protostellar objects (Mardones et al. 1997, henceforth M97); Gregersen et al. 1997; 
Bontemps et al. 1996; Ohashi et al. 1997). Incoherent motions have no directionality and 
appear random, and include thermal motions within molecular gas and nonthermal motions 
of turbulence (i.e., randomly propagating MHD waves; Arons & Max 1975). 

Understanding gas motions is crucial to understanding the overall flow of mass during 
the star formation process, and this is typically attempted by analyzing molecular gas 
line profiles. In practice, however, different types of motion can be difficult to distinguish 
because similar line profiles can be produced when different molecular gas motions are 
projected along the line-of-sight. For example, asymmetrically-blue optically-thick lines 
relative to symmetric optically-thin lines have been interpreted recently as evidence for 
infall motions in several molecular cloud cores (e.g., see Leung & Brown 1977; Zhou et al. 
1993). However, outflow or rotational motions of molecular gas can also produce such a 
signature, especially at velocities that are low relative to the systemic velocity (see Adelson 
& Leung 1988; Cabrit & Bertout 1986). In addition, thermal and nonthermal motions 
each act to broaden line widths, and together can potentially obscure signatures of other 
motions. 

Analysis of the gas motions associated with star formation is further complicated 
by the fact that stars form in groups. For example, young stars are found by number 
more in close proximity within massive, turbulent cloud cores (as in Orion, Perseus, or 
Ophiuchus) than in isolation in less-massive, non-turbulent cores (as in Taurus or Lupus; 
see Lada, Strom, & Myers 1993). Given the flnite resolutions of millimeter telescopes, the 
apparent proximity on the sky of several protostellar objects within a group can confuse the 
relationships between observed molecular gas motions and individual objects, especially in 
more distant regions. 

Millimeter interferometers can minimize the problems of ascribing various types of 
motion to gas within regions where groups of young stellar objects are forming, by providing 
high-resolution observations of judiciously-chosen molecular lines. For example, lines 
observed at high resolution can yield proflles that are less ambiguous in interpretation. 
Furthermore, the high resolution can allow a better determination of the spatial association 
of gas with speciflc objects in a crowded fleld. In addition, high-resolution observations can 
reveal strikingly the abundance contrasts between outflowing and non-outflowing gas, so 
observing species particularly abundant in either of these components can better trace its 
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specific motions. 

In this paper, we present millimeter interferometer observations of the NGC 1333 
IRAS 4 region in the L1430 core of Perseus to probe the individual coherent and incoherent 
motions associated with a small group of protostellar objects. NGC 1333 itself is a region 
forming many protostellar groups, and has a 350 pc distance^ (Herbig & Jones 1983). 

IRAS 4 consists of 4 (or more) very-embedded objects located within a 30" (0.05 
pc) radius, all ~3' (0.3 pc) southeast of SVS 13. IRAS 4 was discovered by Jennings et 
al. (1987) as a single IRAS source, and later resolved by Sandell et al. (1991) into 4A 
and 4B, separated by 30". Although 4A and 4B together were considered one of the first 
"protobinaries" ever detected (see Mathieu 1994), each object is itself multiple. IRAS 
4A was resolved into 4A and 4A', with 2" separation, by Lay, Carlstrom & Hills (1995), 
and IRAS 4B was resolved into two objects, 4B and 4B', with 10" separation, by Looney, 
Mundy & Welch (2000). The 4B object itself also may be a binary or triple system of <0"5 
separation, although this remains uncomfirmed. The 4B' object has been named differently 
by other groups who detected it contemporaneously, i.e., 4BII by Smith et al. (2000), 4BE 
by Sandell & Knee (2000), and 4C by Looney et al. and Choi, Panis & Evans (1999). We 
name the object 10" east of 4B as 4B' because 4C was used by Rodriguez, Anglada & Curiel 
(1999) and Sandell & Knee to name a millimeter object ~40" east-by-northeast of 4A. The 
spectral energy distributions (SEDs) of 4A and 4B are both very red, suggesting that IRAS 
4 consists of several "Class 0" objects, each at the earliest observed stage of protostellar 
evolution (i.e., t ^ 10^ years or less; see Andre, Ward-Thompson & Barsony 1993). 

IRAS 4 has been identified previously as a site of outflow and possibly infall in a 
turbulent region forming clusters. Using single-dish observations of numerous millimeter 
and submillimeter emission lines, Blake et al. (1995, henceforth B95) probed abundance 
depletions and enhancements in the gas associated with 4A or 4B. Low velocity line emission 
(i.e., at AV^ fa 1-2 km s~^ relative to the systemic velocity) revealed a non-outfiowing gas 
component around each source, with abundance depletions by factors of 10-20 (relative 
to dust) due to the adsorption of molecules onto grains at the high densities (n = 2 x 
10® cm~^) and low temperatures (T^ = 20-40 K) of the gas surrounding each object. The 
profiles of some optically-thick lines seen by B95 were asymmetrically-blue, suggesting these 
gas components were infalling. Outflowing gas around each source was also detected from 
line wings seen at higher relative velocities. Qualitative differences in the outflows were 
found between sources, e.g., the 4A outflow extended 4' tip-to-tip in CO and 1' tip-to-tip 



''A distance to NGC 1333 of 220 pc has been suggested by Cernis 1990, but we assume the 350 pc value 
for this paper. 
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in CS, but the 4B outflow extended only <20" in either tracer. Despite these differences, 
however, abundance enhancements (relative to CO) were found in both outflows for several 
molecules including H2CO, CS, and CH3CO, likely due to the return of grain mantles to 
the gas phase from grain-grain collisions within outflow shear zones. 

Further evidence for infall motions at 4A and 4B were obtained from the asymmetrically- 
blue proflles of H2CO 2i2-lii and CS 2-1 by M97, and HCO+ 4-3 and HCO+ 3-2 by 
Gregersen et al. (1997). Mardones (1998, henceforth M98) further mapped IRAS 4 in 
H2CO 2i2^1ii and H2CO 3i2-2ii and found asymmetrically-blue proflles throughout the 
region, but the most extreme self-absorptions were found at the positions of 4A and 4B 
(see also Mardones et al. 2001). Choi, Panis & Evans mapped IRAS 4 in HCO+ 1-0 and 
HCN 1-0, and again found asymmetrically-blue emission proflles at the positions of each 
protostellar object. In addition, small redshifted absorption dips were seen in their HCN 
spectrum toward 4A, arising from infalling gas along the line-of-sight. 

We describe our millimeter interferometer observations in §2. A description of the 
H2CO 3i2^2ii, CS 3-2, and N2H+ 1-0 data is presented in §3. Discussion of the infall and 
outflow motions revealed by H2CO and the rotational and turbulent motions revealed by 
N2H+ follows in §4. In §5, a concluding summary is presented. An Appendix describes the 
continuum data and the relative millimeter SEDs of 4A, 4B, and 4B'. 



2. Observations 

IRAS 4 was observed between July 1997 and March 1998 at the IRAMQ Plateau de 
Bure Interferometer (PdBI), NROQ Nobeyama Milhmeter Array (NMA), and the NRAOQ 
Very Large Array (VLA). Table 1 lists the wavelengths, frequencies, conflgurations, and 
dates of observations made with each array, as well as the synthesized beam FWHMs and 
the 1 a (rms) continuum sensitivities obtained from combining all respective tracks. 
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2.1. Primciry Observations 

2.1.1. PdBI Data 

At the IRAM PdBI, H2CO 3i2-2ii at 225.7 GHz was observed to trace infall or 
outflow motions in the region, following M97 or B95 but using other H2CO transitions 
and at higher resolution. Dense gas can be traced by this line, as its critical density ricr 

7 X 10^ cm^^ (from Green 1991), but H2GO may be depleted in dense gas surrounding 
protostellar objects and enhanced in outflows (sec B95). In addition, N2H+ 1-0 at 93.17 
GHz was observed to trace non-outflowing dense gas in IRAS 4, following M97 but at higher 
resolution. This line is abundant in star-forming regions (e.g., Womack, Ziurys & Wyckoff 
1992) and has also a high critical density, i.e., Ucr ~ 2 x 10^ cm~^. Furthermore, N2H+ is 
observed and predicted to remain relatively abundant in dense cores with n ^ 10^~^ cm~^ 
(Benson, CaseUi, & Myers 1998, Bergin & Langer 1997), but is possibly depleted, rather 
than enhanced, by outflow chemistry (see Bergin, Neufeld & Melnick 1998). Continuum 
emission at A = 3.2 mm and A = 1.3 mm was observed to locate each protostellar object in 
IRAS 4 from associated warm dust. 

The 1 mm and 3 mm bands were observed simultaneously. Receiver 2 was tuned to 
observe H2CO 3i2-2ii over 230 channels in a 20 MHz wide band, yielding a channel spacing 
of 0.10 km s~^ or a velocity resolution of 0.13 km s~^. The remaining 2 correlator units 
were assigned to observe continuum emission at A = 1.3 mm over a 300 MHz wide band. 
Receiver 1 was tuned to observe N2H"'" 1-0 over 230 channels in a 10 MHz wide band, 
allowing all 7 of its hyperflne components to be sampled with a channel spacing of 0.13 km 
s~^, or a velocity resolution of 0.16 km s~^. The other 2 correlator units were combined to 
observe continuum emission at A = 3.2 mm in a 300 MHz wide band. 

Each track with the IRAM PdBI consisted of 3 pointings centered on 4A, 4B, and 
a position halfway in-between, interleaved with observations of 0415+379 (3C 111) or 
0316+413 (3C 84) for phase and flux calibration. (The sources 2230+114 and MWC 349 
were also observed as secondary flux calibrators.) With primary beam FWHMs of ~54" 
and ~22" at 93.17 GHz and 225.7 GHz respectively, this pointing strategy allowed the 
construction of an over-sampled mosaic at A = 3.2 mm and an under-sampled mosaic at 
A = 1.3 mm. Each track was calibrated in Grenoble, France using the IRAM software 
package GLIC. Data from each pointing were combined in the visibility plane and inverted 
simultaneously, and the resulting mosaics cleaned using the IRAM software package 
MAPPING. Further analysis was carried out using the MIRIAD, GLASS, and AIPS 
software packages. 
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2.2. Secondeiry Observations 

2.2.1. NRO NMA Data 

The CS 3-2 line at 146.97 GHz was observed at the NRO NMA in order to trace dense 
gas, as its ricr ~ 1-3 x 10^ cm~^ (Evans 1999). Continuum emission at A = 2.2 mm and A 
= 2.0 mm was observed in order to trace warm dust, improving the SED of each source. 

The NMA FX correlator was tuned to observe CS 3-2 at 146.97 GHz, over 1024 
channels in a 32 MHz wide band, for a channel spacing and velocity resolution of 0.064 km 
s~^. The NMA Ultra- Wide Bandwidth Correlator was tuned over a 1024 MHz wide band 
to observe simultaneously continuum emission at A = 2.2 mm and A = 2.0 mm in the lower 
and upper sidebands respectively. 

Each track with the NRO NMA consisted of 2 pointings centered on 4A and 4B, 
interleaved with observations of 0316+413 (3C 84) for phase and flux cahbration. With 
a primary beam at 146.97 GHz of ~51", this pointing strategy resulted in slightly 
under-sampled mosaics at A = 2.2 mm and A = 2.0 mm. Data from each track were 
calibrated in Nobeyama, Japan, using the software package UVPR0C2. Further analysis 
was carried out using the MIRIAD software packages. For better sensitivity of CS 3-2, 
the NMA data were tapered with a 4" x 4" FWHM Gaussian to obtain a final 6'.'5 x 4'.'4 
FWHM synthesized beam. 

2.2.2. NRAO VLA Data 

Continuum emission at A = 6.9 mm and A = 1.3 cm was observed with the NRAO 
VLA to trace dust, extending the SED of each source across the entire millmeter wavelength 
range. Observations at these wavelengths are also important because they set effective 
hmits on the amounts of free-free emission within the SEDs of young stellar objects at 
shorter millimeter wavelengths, allowing the amounts of thermal emission solely from dust 
to be characterized more accurately (e.g., see Di Francesco et al. 1997). 

Both bands were observed simultaneously by dividing the VLA into two sub-arrays, 
with 13 antennas observing at A = 6.9 mm (43.34 GHz) using Q-band receivers, and 14 
antennas observing at A = 1.3 cm (22.46 GHz) using K-band receivers, each over 100 MHz 
wide bands. 

The NRAO VLA was pointed at 1 position halfway between 4A and 4B to include 4A, 
4B, and 4B' within the ~62" and ~120" primary beam FWHMs of the VLA antennas at 
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A = 6.9 mm and A = 1.3 cm respectively. Observations of IRAS 4 were interleaved with 
observations of 0316+413 (3C 84) for phase calibration approximately every 10 minutes 
and X-band reference pointing checks approximately every 60 minutes. Observations 
of 0134+329 (3C 48) and 1328+307 (3C 286) were made for flux calibration at both 
wavelengths. Data from this track were calibrated in Socorro, U.S.A., and further analyzed 
using the software package AlPS. 

3. Results 

Figure 1 shows maps of IRAS 4 in continuum emission at A = 1.3 mm, A = 2.0 mm, 
and A = 3.2 mm, as well as maps of the integrated intensities of the H2CO 3i2-2ii, CS 
3-2, and N2H+ 1-0 line emission (i.e., the "zeroth-moment" obtained at a given position by 
summing emission from all channels except those containing values < 2 a.) 

3.1. Continuum 

Figures la, lb, and Ic illustrate the relative positions (in projection) of protostellar 
continuum objects in IRAS 4, namely 4A, 4B, and 4B' at 3 different wavelengths. These 
objects were also detected at >4 c at A = 2.2 mm and A = 6.9 mm, but only 4A and 4B 
were detected at A = 1.3 cm. Table 2 lists flux densities measured from our data (also see 
Appendix.) 

Separated from 4A by only ~1", 4A' is not detected distinctly from 4A in our ~2-4" 
resolution maps, but the continuum emission around 4A is extended towards its position. 
For this paper, we limit discussion of 4A and 4A' as an unresolved pair rather than consider 
them individually, i.e., 4A = (4A + 4A'). Similarly, we refer to the 4B object as 4B, though 
it may also be multiple. Furthermore, the nearby source 4C was not observed as it lay 
outside the fields-of-view of our PdBI and NMA mosaics, and our VLA images did not have 
sufficient sensitivity to detect it (based on its flux densities at A = 3.6 cm and A = 6.0 cm; 
see Rodriguez et al.) 

3.2. H2CO 3i2-2ii and CS 3-2 

Figures Id and le show that H2CO and CS are distributed similarly around 4A, 4B, 
or 4B', but remarkable differences are found between objects. Near 4A, an extended lobe is 
located to the south-to-southeast and 1 or 2 compact spots are located to the north in both 
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H2CO and CS. The southern lobe contains 4 local maxima in the H2CO data, but only 2-3 
maxima in the lower-resolution CS data. Near 4B, 2 single-peaked, compact lobes are seen 
to the north and south, with similar brightnesses in H2CO but very different brightnesses 
in CS. In addition, three less prominent peak;s to the east, northwest and southeast of 4B 
are detected in H2CO but not in CS. Finally, no H2CO or CS emission is detected near 4B'. 
In our H2CO map, this source is located at the mosaic edge, making it unclear if the lack of 
detected H2CO near 4B' is real or due to reduced sensitivity. In our CS map, however, 4B' 
is positioned well within the fields-of-view, but no CS is detected near the object. 

Figure 2 shows grids of H2CO 3i2-2ii spectra from positions centered on 4A and 4B, 
spaced regularly in 2" increments of R.A. or Dec, i.e., ~1 FWHM of the synthesized beam. 
The velocity range of each panel is centered at -1-6.96 km s~^, the average central velocity 
determined from N2H+ 1-0 data of 4A and 4B (see §4.3.2 below), which we consider here 
to be the systemic velocity of the IRAS 4 group. 

Figure 2a reveals that at the 4A position (i.e,. the central grid panel), an "inverse 
P-Cygni" hue profile of blue-shifted emission and red-shifted absorption is clearly detected. 
The absorption feature has a central velocity of 7.6 km and a width of ~0.9 km s~^, 
and is not prominent at any other grid location. The maximum depth of the absorption 
toward 4A is ~9 K, or 1.3 Jy beam~^, approximately equal to the peak intensity of A = 
1.3 mm continuum emission detected at the same position (see Table 2). With single-dish 
data, M98 detected asymmetrically-blue, single-peaked H2CO emission with no absorption 
at the 4A position. North of 4A, the H2CO spectra are double-peaked, with a dip at ~7.7 
km s~^, and a weak red wing is seen. South of 4A, the H2CO spectra are single-peaked, but 
a strong blue wing is seen. The emission spot due north of 4A in Figure Id almost entirely 
consists of this red-shifted wing emission. (The apparent H2CO emission spot northwest 
of 4A is incidental map-edge noise.) The extended lobe seen south-to-southeast of 4A in 
Figure Id almost entirely consists of blue-shifted wing emission. 

Figure 2b reveals an inverse P-Cygni hue profile at the 4B position (at the center of 
this grid, defined similarly to Fig. 2a). Here, the absorption has a central velocity of 7.9 km 
s"^ and a width of ~1.0 km s~^. The maximum depth of the absorption feature toward 4B 
is ~4 K, or 0.57 Jy beam~^, again approximately equal to the peak intensity of continuum 
emission at A = 1.3 mm. In addition, absorption is not prominent at positions offset from 
the continuum object. M98 also only detected asymmetrically-blue, single-peaked H2CO 
emission with no absorption below the continuum at the 4B position. Furthermore, spectra 
to the north of 4B are double-peaked and a red wing is seen, like that seen north of 4A and 
also with a self- absorption at ~7.7 km s~^. South of 4B, the hue is single-peaked and a blue 
wing is found, similar to the line south of 4A. The bright peaks north or south of 4B seen 
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in Figure Id predominantly consist of either red- or blue-shifted wing emission respectively, 
and the secondary features seen northwest and soTithcast of 4B consist of blue-shifted and 
red-shifted emission respectively. (The third emission feature seen in Figure Id near 4B is 
off the eastern edge of the grid in Figure 2b, and has both blue- and red-shifted emission.) 

Line profiles similar to those found for H2CO 3i2-2ii in Figure 2 are also seen for the 
CS 3-2 line at similar positions, including an absorption feature at 7.8-8.0 km s~^ in the 
CS spectrum at 4A, and a possible analogue at 4B. These data arc not shown since they 
corroborate the H2CO data, but are noisier and have lower spatial resolution. 

3.3. N2H+ 1-0 

Figure If shows the integrated intensity (zeroth-moment) of N2H+ 1-0 emission around 
the IRAS 4 objects. These distributions differ significantly from those of H2CO and CS 
(Figures Id and le) around each object. However, the spatial distribution of the integrated 
intensity also differs from object to object. Near 4A, the N2H+ emission is elongated 
northeast-southwest, extending somewhat symmetrically 10-15" around a bisecting channel 
at position angle (P-A.) fa 45° containing 4A and 4A'. Near 4B, the N2H+ emission is more 
compact than that near 4A, although it is slightly elongated north-south. In addition, 
N2H+ emission is coincident with the 4B object, unhke at 4A (i.e., no bisecting channel is 
seen in the zeroth-moment map.) Near 4B', N2H+ is not detected, although its position 
here is well within the field-of-view, unlike in the H2CO mosaic. 

Figure 3 shows grids of spectra of N2H+ 101-012, at the same positions at and around 
4 A and 4B shown in Figure 2. (For clarity, only the "isolated" hyperfine component of 
N2H+ 1-0 is shown.) Since the 2" spacing of panels in Figure 3 corresponds approximately 
to half-beam spacing at the FWHM of the synthesized beam at 93.2 GHz, neighboring 
panels are not independent, as in Figure 2. (In addition, the velocity range of these grid 
panels differs from that of Figure 2, but remains centered at the same systemic velocity.) 

Figure 3a reveals (in the central panel) another inverse P-Cygni line profile at the 4A 
position, with absorption of central velocity and width similar to those of H2CO at the same 
position, i.e., 7.9 km s~^ and ~1 km s^^ respectively. (Absorption seen here in adjacent 
panels is still coincident with the continuum object due to the half-beam, non-independent 
grid spacing.) Unlike the H2CO absorption, the peak depth is only ~1.2 K, or ~50% the 
peak continuum intensity of 4A at A = 3.2 mm. M98 did not detect absorption of N2H"'" 
1-0 in their single-dish data at the 4A position. At positions northwest and southeast of 
4A, the line is single-peaked in emission only with small shifts in central velocity to the red 
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and blue respectively. No systematic line width variations or wing emission is seen. The 
bisecting channel at a P.A. 45° seen in Figure If is noticeable here from the relatively 
dim emission found in panels on the northeast-southwest diagonal. 

Figure 3b shows that an inverse P-Cygni line profile is not seen at the position of the 
4B object, rather a fairly symmetric line is found peaking at 7.1 km s~^, near the systemic 
velocity. (Absorption may be present at ~8 km s~^ in the grid panel 2" west of the central 
position, but this is not suggested in any other adjacent, non-independent grid panel.) 
Unlike near 4A, there does not appear to be substantial variations in central velocity over 
the grid, but the line is dimmer and broader south of 4B and brighter and narrower north of 
4B. Like near 4A, there is no evidence for wing emission. M98 detected a similarly-shaped 
single-peaked line profile at the 4B position. 

4. Discussion 

4.1. Distinct Tracers of Specific Motions 

The spectral and spatial distributions of the lines described in §3 trace quite distinctly 
specific types of motion within IRAS 4. Infall is traced from the inverse P-Cygni profiles 
detected in H2CO toward both 4A and 4B (and also in CS and N2H+ toward 4A alone.) 
Outflows from 4A and 4B are traced distinctly from infall, at positions offset from the 
protostellar objects, from H2CO and CS lines with profiles of only bright line wings and 
no line cores. Rotation of non-outflowing dense gas around 4A or 4B is traced distinctly 
from outflow from spatial variations of the central velocity of N2H+, with proflles of only 
line cores and no line wings. Furthermore, variations of turbulent motions within the 
non-outfiowing dense gas about 4A and 4B are traced by spatial variations of the fine width 
of N2H+ in the same profiles. 

The observed lines distinctly trace specific motions partly because of the differing 
gas-grain chemistry of the molecular species. B95 suggested H2CO and CS were depleted 
in the non-outfiowing gas associated with 4 A and 4B, due to significant amounts of these 
molecules adsorbing onto dust grains. Detailed models by Bergin & Langer (1997) of 
relative abundances in centrally-condensed cores of increasing density support this scenario, 
showing H2CO or CS may deplete by >2 orders of magnitude at densities of n ~ 10^ cm~^ 
onto grains with H2O or CO mantles. These models, in turn, also suggest that N2H"'" 
remains at a constant abundance to densities of ~10^ cm^^, due to the low binding energy 
of N2, its precursor molecule, to grain mantles. Therefore, we expect that N2H^ traces 
non-outfiowing dense gas better than H2CO or CS in a suitably dense protostellar core. 
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In contrast, B95 suggested the enhanced abundances of H2CO and CS found in the hne 
wings were due to significant amounts of these molecules returning to the gas phase from 
the disruption of grain mantles in the outflows. Moreover, the single-dish data of M97 
show the N2H^ 1-0 components at 4A and 4B have Gaussian profiles, with no emission 
at the wing velocities of H2CO or CS. Models of molecular abundance variations from 
post-shock chemistry by Bergin et al. support this contrast, since the increased abundance 
in outflows of H2O, returned to the gas phase from grain mantles, appears to rapidly deplete 
N2H"'". Therefore, we expect H2CO or CS traces better outflowing gas than N2H''". Such 
molecular abundance contrasts are not unique to IRAS 4 and have been observed similarly 
elsewhere, e.g., the single-dish observations of the LI 157 core and outflow by Bachiller & 
Perez-Gutierrez (1997). 

Another important reason our data distinctly trace specific motions is their high spatial 
resolution. For example, interferometers allow detection of the inverse P-Cygni profiles 
because the brightnesses of the continuum objects in the relatively small beams can exceed 
the brightnesses of line emission from intervening, infalling gas along the line-of-sight, 
leading to red-shifted absorption. (Blue-shifted emission in the inverse P-Cygni profiles 
arises from concentrated gas on the far sides of the protostellar objects.) In addition, the 
spatial filtering properties of the interferometer can further distinguish specific motions 
better than single-dish observations, especially at low relative velocities, by selectively 
sampling the kinematics of compact structures, and resolving out confusing, extended 
emission. For example, outflowing gas is detected more distinctly in interferometer 
observations of H2CO or CS, beyond abundance differences, because zones of outflowing 
dense gas are expected to be more compact, and hence are more easily seen by an 
interferometer, than the outer, less-dense layers of a circumstellar envelope that may 
still contain abundant H2CO or CS. On the other hand, interferometers can still detect 
non-outflowing gas in such an envelope using N2H+ because, unlike H2CO or CS, it remains 
relatively abundant in the denser and more-compact inner regions that are less resolved 
out than outer regions. (As described below, N2H+ may be somewhat depleted relative to 
"canonical" values in IRAS 4.) 



4.2. H2CO Observations 

4.2.1. Infall 

Detections of red-shifted absorption indicative of infall motions have been made 
primarily toward high-mass young stellar objects, where the absorption profiles were seen 
against the continuum emission of bright HII regions in interferometer beams. For example. 
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red-shifted absorption was seen toward GlO.6-0.4 in NH3 (1,1) and (3,3) by Keto, Ho 
& Haschick (1987, 1988; with VIA beams of 3"0 FWHM and 0'.'3 FWHM respectively). 
Similar profiles were also detected toward W51:e2 in NH3 (1,1) and (2,2) by Ho & Young 
(1996; with VLA beams of 2'/6-3'.'8 FWHM), and in these transitions plus NH3 (3,3) by 
Zhang & Ho (1997; with VLA beams of 0'.'2-l'/0 FWHM), and in CS 3-2 by Zhang, Ho & 
Ohashi (1998; with an NMA I'.'l FWHM beam). Furthermore, an inverse P-Cygni profile 
was noted toward W49A North: G in HCO+ 1-0 by Welch et al. (1987; with a BIMA ~7" 
FWHM beam). 

Toward low-mass young stellar objects, detections of red-shifted absorption against 
continuum emission have been very rare, the only case being the weak inverse P-Cygni 
profile detected in the HCN 1-0 line toward NGC 1333 IRAS 4A by Choi, Panis & Evans 
with a BIMA ~8-12" FWHM beam. (Red-shifted absorption is also suggested in their 
Figure If in the HCN line toward IRAS 4B, and in the HCO+ 1-0 line toward both objects, 
but they describe only the HCN line toward 4A as an inverse P-Cygni profile. )0 

Figure 4 shows the H2CO 3i2-2ii line at positions of maximum continuum intensity 
from 4A and 4B (enlarged from the central panels of Figures 2a and 3a.) These data provide 
the least ambiguous evidence for infall motions toward low-mass young stellar objects. 
First, the detection of the absorption here is very strong, given the low channel-to-channel 
noise of these data. Second, the excitation requirements of the lower energy level of the 
H2CO 3i2-2ii transition (i.e., -^211/^ = 7.5 K) are higher than those of HCN 1-0, suggesting 
the absorptions are more likely occurring in warmer gas associated with the young objects, 
rather than in a layer of cooler gas at slightly higher velocity situated along the line-of-sight 
(and possibly not associated with the young objects.) In addition, the scales probed by 
these data are smaller than other investigations. At a distance of 350 pc, the linear beam 
size (FWHM) of the H2CO data at IRAS 4 is ~0.003 pc, a factor of ~4-6 less than the 
scales probed by Choi et al., and a factor of ~2 less than those probed with inverse P-Cygni 
profiles at the more-distant, higher-mass objects noted above. 

To probe conditions in the absorbing gas, the inverse P-Cygni profiles toward 4A 
and 4B were modeled with an enhanced version of the "two-layer" code used by Myers 
et al. (1996) to reproduce successfully the appearance of spectral lines from infalling gas 
toward several protostellar cores (see also Williams & Myers 2000, and Lee, Myers & 
Tafalla 2001) In this version, two layers, located along the line-of-sight, each have peak 
optical depth Tq, velocity dispersion a, and approach speed Vi„ toward a continuum source 



^An inverse P-Cygni profile of the ^'^CO (2-1) line toward 4A has been also recently detected by Webster 
& Welch (2000) with a BIMA 6'.'7 FWHM beam. 



-14- 



located in-between. The continuum source is assumed to be optically thick with a Planck 

temperature J,., where is related to the blackbody temperature at frequency by Jc = 
To/(exp(To/Tc)-l), with = hu/k, h = Planck's constant, and k — Boltzmann's constant. 
The continuum source also fills a fraction $ of the telescope beam area. The "front" layer 
between source and observer has a Planck excitation temperature Jf while the "rear" layer 
behind the source has a Planck excitation temperature J^. The rear layer is illuminated by 
the cosmic background radiation of Planck temperature J^. 

For this system, the observed line brightness temperature at velocity V can be written 

as 

ATb = {Jf - Jcr)(l - exp(-r/)) + (1 - *)(J, - Jb){l - exp(-T, - t/)), (1) 

where 

Je. = $Jc+(l-$)J., (2) 

and, 

Tf = Toexp{-{V - - VLSR)V2a'), (3a) 

Tr = T„exp(-(y + Vin - Vi^SR? jlo'') . (36) 

In this model, the line profile can be understood very simply from Equation 1 as the 
linear sum of a broad positive profile centered on Vlsr (i-e., the right term of RHS) and a 
narrow negative profile centered on VLSi^+P^n (i.e., the left term of RHS). If the systematic 
motion of the layers is inward (i.e., Vin > 0), Equation 1 implies that the resulting profile 
shows red absorption or "infall asymmetry." If the front layer is cooler than the continuum 
source and rear layer by a relatively small temperature differences Jf-Jcr, the asymmetric 
profile shows "self-absorption" but remains positive at all velocities. As this temperature 
difference increases, the profile becomes more asymmetric and its "dip" goes below the 
spectral baseline. 

Equations 1-3 are essentially the same as Equations 1 and 2 of Myers et al. (1996), 
except that the emission from an optically thick central source is now included. This 
model is intended to give the simplest possible description of a contracting or expanding 
system which can reproduce the observed line profiles. Therefore, it gives only the crudest 
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picture of the optical depth and excitation temperature of the absorbing and emitting gas. 
However, the model is still useful for comparing the characteristic velocity Vin from one line 
profile to the next, since the ratio Vin/cr is most important to setting the line profile for any 
model of a contracting system (Leung & Brown 1977). 

In addition to the observed spectra. Figure 4 shows model inverse P-Cygni profiles 

that successfully reproduce those observed. The respective model parameters arc listed in 
Table 3. Infall speeds of Vin = 0.68 km s^^ and 0.47 km s~^ are obtained toward 4A and 4B 
respectively. (The blue wing in the 4A spectrum is likely from outflow and so it is not fit by 
the model.) Using essentially the same two-layer model (but without a continuum object), 
M98 found a V^n of only 0.11 km s~^ for the IRAS 4 region from the same H2CO line. 
The difference between these results is hkely because the spectrum modeled by M98 was 
obtained with a resolution coarser by a factor of ~6 than those shown here, and consisted 
of data averaged over the FHWM of their single-dish integrated intensity map, i.e., over a 
projected distance from the source of ~0.04 pc. The resulting emission spectrum included 
data from positions offset from 4A and 4B with gas having possibly lower apparent infall 
velocities due to projection, whereas the inverse P-Cygni profile traces material along the 
line-of-sight to 4A or 4B with minimal projection effects. 

The inverse P-Cygni profiles alone cannot demonstrate whether the infalling material 
consists of dense gas individually surrounding 4A or 4B, or from a common layer surrounding 
both objects, since infalling material at projected locations between 4A and 4B is not 
traced. However, the difference between Vm toward 4A and 4B is significant in that the 
range of Vin able to reproduce well each spectrum is only a few hundredths of km s~^, if 
all other parameters are kept fixed. Different values of Vin toward 4A and 4B suggest the 
absorptions arise in dense gas residing around each source rather than in a layer common to 
both objects, where presumably more similar velocities toward 4A and 4B would be found. 

From simple gravitational arguments, the line-of-sight speeds implied by the absorptions 
allow crude estimates of the masses of the 4A and 4B objects, assuming the infall speeds 
Vin arise from the velocity gain of gas free-falling from rest at r = 00 to r = rj„, or 

M ^ (V.^.„)/2G ^ ,.2[^^J {0-)m,. (4) 

To obtain (an upper limit of) rj„, the mass obtained from Equation 4 can be equated with 
that within a sphere of radius rj„ of mean density n, i.e., M = {4/3)7irf^^nfimH , with fi equal 
to the mean molecular weight (2.32), and ttih equal to the mass of hydrogen. We assume 
n = 2.3 X 10^ cm~^, the critical density of the H2CO 2ii-lio line, a sufficient density for 
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substantial population of the 2ii energy level of H2CO, determined at = 30 K (Green 
1991) since B95 estimated a kinetic temperature of non-outfiowing gas around 4A of 20-40 
K. Accordingly, the values of Vin obtained for 4A and 4B yield values of rj„ of 0.0099 pc 
and 0.0068 pc respectively, or 5.8" and 4.0" respectively at the assumed distance of IRAS 
4. These radii are smaller by factors of 1.4-2.2 than those obtained using integrated N2H"'" 
1-0 emission (see Figure If) to trace the projected extents of dense gas around 4A and 4B, 
which is consistent given the lower critical density of that transition.^ Using Equation 4, 
the masses interior to the respective values of rj„ are 0.53 M0 and 0.17 M© for 4A and 4B, 
reasonable values given the moderate luminosities of these objects and their location within 
a region of low-mass star formation. 

Mass accretion rates. Mm, can be estimated toward 4A and 4B using the simple 
expression: Mj„ = inrfj^fimnnVin, or 

- ^ ^°"fe)'( 2.3xl"o.cn.-j (l^) • 

assuming spherical symmetry of the infall. With the values of Vi„, n, and rj„ described 
above, we estimate = 1.1 x 10""^ Mq yi'^ for 4A, and 3.7 x 10"^ M© yr'^ for 4B, 
using Equation 5. Assuming these rates are constant, accretion timescales tacc can be then 
estimated by dividing the masses found above with the accretion rates. Such a calculation 
yields tacc = 4.7 x 10'^ yr and 4.6 x 10^ yr for 4A and 4B respectively, remarkably similar 
and well within the tacc < 10^ surmised for Class objects. 

If the isothermal sound speed of this gas can be approximated as K-ms = {kT^/ fimn)^ 
and Tfc ^ 30 K, the values obtained here for V^„ suggest the infall motions toward 4A or 
4B are supersonic, by factors of 2.1 and 1.4 respectively. Given that the magnetic field 
strength is poorly known in IRAS 4, however, it is not clear if the infall motions are also 
sub- or super-Alfvenic. From a virial theorem argument, Akeson et al. (1996) estimated 
large magnetic field strengths associated with 4A of 5-8 mG, suggesting the infall flows are 
very sub-Alfvenic. However, the summary of actual magnetic field measurements (of other 
various molecular clouds) by Crutcher (1999) suggests these values may be too large by 
factors of ~10-100. If the infall motions are occurring exactly at the Alfven speed (i.e., Va 
= B / {ATCfinmfjY^'^ with n = 2.3 x 10^ cm^^, the required magnetic field strengths near 4A 
and 4B would be 720 /iG and 500 fxG respectively. 



^Non-outflowing H2CO in emission, i.e., over the velocities traced by N2H+, is difficult to see around 4A 
and 4B, possibly due to low abundances or the relative brightness of the outflows. 
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The two-layer models used to compute the line profiles seen in Figure 4 neglect 
the radial variations in density, temperature, and velocity expected within protostellar 
envelopes. Therefore, these results are too simple to be definitive, but provide motivation for 
more observations and models of inward motions. It will be useful to make more-detailed, 
Monte-Carlo models of radiative transfer through envelopes, of the H2CO line from these 
and other Class objects. 

4.2.2. Outflow 
4.2.2.1 Outflow Morphologies and Origins 

Figure 5 shows channel maps of H2CO 3i2-2ii emission near 4A, revealing the 
underlying structure of the blue 4A outflow lobe better than the integrated intensity 
map (Figure Id). In channel maps, the blue lobe consists of two prominent features 
south- by-southwest and southeast of 4A. Increasing in relative velocity, emission from the 
first feature is first seen at 6.5 km s~^, and it peaks at increasing angular separations from 
4A in a Hubble Law manner reminiscent of other protostellar outflows (see Bachiller 1996). 
At 2.0 km s^^, however, the angular separations of the peaks of this feature shift to the 
southeast with increasing velocity out to the band edge, -10 km s^^. The second feature is 
first seen at 6.8 km s^^ and it moves south and west with increasing relative velocity in a 
manner less consistent with a Hubble Law. Together, these features resemble cavity walls of 
an outfiow seen at increasing relative velocities. In contrast to Figure Id, the channel maps 
suggest the 4A blue lobe has an overall P. A. of ~160°, and its north-south appearance in 
Figure Id is due to the relative brightness of the western cavity wall. Figure 5 also shows 
dense gas outfiowing from 4A', the binary companion to 4A, at low relative velocities along 
the south-by-southwest feature, but dense gas appears to be outfiowing instead from 4A at 
higher relative velocities. 

Our H2CO data reveal a systematic decrease in the P.A. of the 4A outfiow axis with 
scale that could be due to precession of the 4A outfiow axis (i.e., a systematic change in 
its P.A. with time.) Single-dish CS data from B95 and Lefioch et al. (1998) show the blue 
lobe oriented southwest at separations >35" from 4A but oriented south at separations 
15-35" from 4A. Our data show this trend continues to smaller separations, since the blue 
lobe is oriented to the southeast at <15" offsets from 4A. Such precession may be due to 
the binarity of 4A, if the gravitational action of 4A' on 4A is strong enough (or vice- versa.) 
The inherent periodicity of precession, however, suggests outfiow gas will be distributed 
periodically in P.A. at even larger scales, and indeed such evidence may exist in recent CO 
maps of NGC 1333. For example, a CO 3-2 JCMT map of Knee & Sandell (2000) shows 
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what may be the northern red lobe of the 4A outflow curving north of IRAS 4, from a P.A. 
of 45° to a P.A. of 0° over 3-4'. In addition, the CO 1-0 BIMA mosaic of NGC 1333 by 
Plambeck & Engargiola (1999) faintly shows possible emission from the 4A outflow at large 
scales, curving symmetrically over several arcminutes from 4A in an S-shape. However, 
maps that are more sensitive than these are required to determine conclusively if this 
distant emission is actually related to the 4A outflow. (As this paper was being revised after 
submission, a periodicity in the P.A. of the 4A outflow was seen by Choi (2001). However, 
the amphtude suggested here is larger by a factor of ~10.) 

The difference in integrated intensity between sides of the 4A blue lobe could be 
understood if the outflow is precessing, as dense gas surrounding 4A may be cleared out 
in some directions more than others by a precessing outflow. For example, the emission 
from the western wall of the outflow cavity shown in Figure 5 may extend to higher relative 
velocities than the eastern wall if the outflow has encountered less mass to the west, a 
direction along which gas may have been most recently cleared out by the precessing 
outflow. This scenario may also explain why the terminus of the outflow cavity appears to 
shift to the west with increasing relative velocity. 

Figure 6 shows channel maps of the dense gas features near 4B, revealing the underlying 
structure of its outflow. Here, all features are found at increasingly large positional offsets 
from 4B with increasing relative velocity in a Hubble Law manner, except the feature 
east of 4B which moves in the opposite sense. Correspondingly, this latter feature may be 
emission from dense gas outflowing from 4B'. (Our mosaic does not cover area east of 4B' 
where other possible outflow features from 4B' could reside.) 

The red features to the north and southeast and the blue features to the south and 
northwest seen in Figure 6 together suggest an X-shape, possibly indicative of outflow cavity 
walls similar to those surmised for the 4A blue lobe. In this case, the 4B dense gas outflow 
may consist of two lobes, one situated at a north- by-northwest P.A. (i.e., ~345°) with an 
opening angle of ~45°, and another at a south- by-southeast P.A. (i.e., ~160°) with an 
opening angle of ~30°. The peculiar appearance of both red and blue outflow emission on 
the same sides of 4B then may be due to the 4B outflow being inchned from the plane-of-sky 
at an angle less than, e.g., the ~30° opening angle of the narrower lobe. An alternative is 
that the brighter and dimmer pairs of symmetric red and blue dense gas outflow features 
are each independent outflows, with a projected difference in the orientation of their axes 
of ~30-45°, driven by unresolved objects within 4B. As with 4A, the H2CO emission may 
still trace cavity walls in this case, but here they may be unresolved. A close companion to 
4B is suggested by the sub-arcsecond resolution continuum data of Lay et al. and Looney 
et al. but this has not been yet confirmed. 
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4-2.2.2 Physical Properties of Outflows 



Physical properties of the outflowing dense gas can be estimated from the H2CO data, 
with some assumptions. To estimate the mass of outflowing gas, local thermodynamic 
equilibrium (LTE) was assumed and line emission was assumed to be optically thin, 
although these assumptions are likely not vahd within dense gas outflows. To obtain 
de-projected relative gas velocities, both outflows were assumed to originate at the Vlsr 
of 4A or 4B (see §4.3.2) and are inclined 10° from the plane-of-sky due to the low spatial 
coincidence of their red and blue lobes. In addition, excitation temperature and H2CO 
abundance values of T^.^. = 50 K and X(H2C0) = 1.8 x 10~^ or 4.3 x 10~^ were assumed 
for 4A or 4B respectively, as derived by B95 from wing emission. 

Re-binning the original data cube to channels of width AV — 0.31 km s~^ 
(approximately the thermal broadening FWHM of H2CO at 50 K, to obtain independent 
velocity channels I), the total gas mass of a given outflow feature k, Mk, can be estimated 
as 



where and fi^/ arc the average flux density and deconvolved area of feature k in channel 
/ obtained from two-dimensional Gaussian fits. (The projected distance to the nearest 
continuum object was also obtained from Gaussian fitting.) About 6 features were identified 
by eye within each of the 4A and 4B outfiows, based on their >4(j peak intensities and 
persistence in several adjacent channels. In addition, Aji is the Einstein A coefficient of the 
3i2-2ii line, u is the line frequency, D is the assumed distance to IRAS 4, ZiT^.^) is the 
partition function of H2CO at T^xi and Qj is the degeneracy of the 3i2 level. 

Masses {Mout) were obtained by summing Equation 6 over all features. Momenta {Pout) 
or kinetic energies (Eout) were obtained by summing the products of mass and relative 
de-projected velocity per channel, or the products of half the mass and the square of the 
relative de-projected velocity per channel respectively over all features. A characteristic 
timcscalc t was determined as the time for mass of a given feature at a given de-projected 
relative velocity to have travelled to a de-projcctcd distance. (Given the extended nature 
of mass in each channel, t should not be mistaken for the actual time that all mass in the 
channel took to travel from the continuum source.) Dividing the momenta and kinetic 
energies per channel by these timescales per channel and summing over channels and 
features yielded total characteristic momentum fiuxes {FH2C0, forces) and mechanical 
luminosities {Lj^^chi or powers) of the outfiows. 
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The total derived masses for the 4A and 4B dense gas outflows are 6.9 x 10~^ Mq and 
1.5 X 10~^ Mq respectively. Table 4 lists these and the other physical properties derived 
for each dense gas outflow. The two brightest features associated with 4A or 4B contribute 
respectively 45% and 43% of Mout- Note that values for 4A listed in Table 4 include only 
the blue lobe of the 4A dense gas outflow, so total properties of the entire 4A dense gas 
outflow may be larger than those listed by a factor of ~2, assuming the flow is symmetric, 
e.g., for a total Lmech of 0.25 yr~^. 

Quantities in Table 4 for the 4A southern blue lobe alone exceed those of the all lobes 
of the 4B dense gas outflow by factors of 5-17. These differences may be signiflcant if both 
outflows have similar low inclination angles, but this is a very uncertain assumption. For 
example, the mechanical luminosity, L^^ch, of the 4 A blue lobe would decrease by a factor 
of ~90 if the inclination of the flow were 45° instead of 10°, due to the intrinsically large 
dependence on tan i. Differences between the 4 A and 4B dense gas outffow characteristics 
could be expected from differences in the photon luminosities of the protostellar sources. 
Millimeter wavelength luminosities of 4A and 4B, obtained by integrating continuum flux 
densities in Table 2 and those at A = 450 fim and A = 850 fim by Sandcll & Knee, are 0.24 
Lq and 0.12 Lq respectively, only a factor of ~2 different. Other factors, therefore, may 
also be relevant. 

Our values of Mout, Pout, and Eout for the 4A blue lobe are factors of ~60-80 less 
than values derived for the same lobe by Knee & Sandell, assuming their inclination angle 
(45°) and distance (220 pc). These differences are likely because the CO 3-2 data of Knee 
& Sandell (from the JCMT with a 14" FWHM beam) have greater sensitivity to large 
quantities of lower density gas (i.e., ncr[CO 3-2] « 10^ cm~^), suggesting dense gas (as 
seen by the H2CO line) is relatively insignficant in terms of the total mass, momentum, 
and energy of the molecular outflow, comprising <2% of these quantities. However, the 
momentum flux and mechanical luminosities we derive are only factors of ~10 less than 
those derived by Knee & Sandell for the 4A blue lobe, indicating the dense gas has greater 
significance in the total force and power of the 4A outflow. 

4.3. N2H+ Observations 

As discussed above, the N2H+ molecular ion can be relatively abundant in non- 
outflowing dense gas around protostellar objects, and the optical depth of its transitions 
can be low in all but the densest cases (e.g., L1544; see Williams et al. 1999). In addition, 
the N2H+ 1-0 transition has 7 hyperflne components from which the total brightness 
temperature {T^°^), the central line velocity (Vlsr), and the line width {AV) can be 
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estimated to good accuracy via simultaneous hyperfine structure fitting. These quantities 
can be used to estimate the mass of non-outflowing dense gas around 4A and 4B, and probe 
its rotational and incoherent motions. 

Figure 7 shows the values of T^*°*, Vlsr, and AV obtained from fitting the hyperfine 
structure of the N2H"'" line around 4A and 4B, using the HFS routine of CLASS (see 
Caselli, Myers & Thaddeus 1995.) Values of AV shown were obtained by deconvolving the 
measured with the 0.16 km s~^ velocity resolution of the PdBI data. Typical errors 
from HFS fitting per pixel for T^*"*, Vlsr, and AV were 0.94 K, 0.015 km s-\ and 0.032 km 
s~^ respectively. The quantities shown in Figure 7 are restricted to those recovered from 
pixels containing data where the peak intensity of the brightest hyperfine component (i.e., 
N2H"'" 123-012) was >5x the rms of the respective spectral baseline. With this criterion, 
regions of dense gas surrounding each object are defined, similar in extent to those obtained 
by the lowest contour of Figure If from an integrated intensity criterion. 

In Table 5, the mean, rms, minimum, and maximum values of T^*°*, Vlsr, and AV 
in the non-outfiowing dense gas surrounding 4A and 4B are summarized. To obtain these 
values, the fields shown in Figure 7 were re-binned into images with pixels 4'.'0 x 3'.'5 in size, 
about the area of the beam, so beam-averaged quantities at spatially independent locations 
could be examined. (Pixels not meeting the criterion for inclusion in Figure 7 were also not 
included in the re-binning.) Around 4B, the mean value of T^*°* is ~60% larger than it is 
around 4A, but the mean values of Vlsr and AV are similar in both cases. However, the 
rms of Vlsr for the non-outfiowing dense gas around 4A is ~4x larger than that for the 
same gas around 4B, a point discussed further in §4.3.2 below. 

4-3.1. Column Densities and Masses 

With quantities obtained from hyperfine structure fitting, the column density of N2H+, 
A^(N2H+), in a given pixel around 4A and 4B can be estimated. Without unique estimates 
of Tea; and Ttot, however, we assume Tex = 10 K >> hB/3k, T^g = 2.73 K.0. In addition, we 
assume the line is optically thin (i.e., T^*°* ~ TexTtot), which is supported by the symmetric 
appearance of the line over the regions shown in Figure 7. Following similar estimates by 
Womack et al., but with these assumptions. 



^"^Tbg is much larger at the positions coincident with the peak continuum intensity, but since this is only 
significant for relatively few pixels, that difference is ignored here. 
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A^(N2H+) 



10^ 3h hiTl"^) e^-^/*^^- 



(7) 



1.06( J + 1) SttV hB (1 - e-WfcTe.) 



where J is the upper rotational level (1), is the dipole moment of the N2H+ molecular ion 

(3.40 Debyc; Green, Montgomery & Thaddeus 1974), B is the rotational constant of N2H"'" 
(46.586702 GHz; Caselli et al.), Ej is the upper level energy (= hv where z/=93. 1734035 
GHz; Caselli et al.), and T^*°*, and come from HFS fits. (Note that units for AV are 
km s~-^.) 

From Equation 7, mean values of iV(N2H+) in the non-outfiowing dense gas around 
4A and 4B of 2.4 x 10^^ cm~^ and 3.9 x 10^^ cm~^ are found respectively. Table 5 also 
lists the rms, minimum, and maximum values for A^(N2H"'") obtained for this gas around 
each protostellar object. (Decreasing to 5 K reduces the values of A^(N2H+) in Table 5 
by only factors of ~5%, but increasing T^x to 20 K increases them by factors of ~40-50%.) 
The mean column densities of non-outflowing dense gas around both objects are similar, 
although that of the gas around 4B is ~60% larger than that around 4A. These values 
are 2.5-5 x larger than the A'"(N2H''') of 0.95 x 10^^ cm~^ obtained from the mean N2H+ 
spectrum of IRAS 4 from single-dish data of M98, under the same assumptions. Maximum 
column densities listed in Table 5 are 5-8 x larger than this latter value. (The iV(N2H"'") 
quantities listed in Table 5 are actually lower limits, given the assumption of Ttot 0, but 
the quantities should not vary much from those listed if Ttot < 1 as expected.) 

With A^(N2H+), the masses of non-outflowing dense gas traced by N2H+, Myas, can be 
estimated by assuming the emission originates from gas local to 4A or 4B and not from gas 
distributed along the line-of-sight. With this assumption. 



where jj,, mn, and D are the same as for Equation 6, but fi^ is the solid angle subtended 
by the dense gas traced by N2H+ around each object, and X(N2H+) is the fractional 
abundance of N2H+. 

For X(N2H"'"), we assume for 4A and 4B the value of 6.0 x 10^^^ derived by B95 
at 4A by comparing the observed column density of N2H+ with that suggested for H2 by 
continuum observations of dust. This value is relatively low compared to mean values found 
in other cloud cores; Womack et al. and Benson et al. found mean values of X(N2H+) of 4 
X lO"^'^ and 7 x 10~^° respectively from various samples (see also Ungerechts et al. 1997). 
However, the sample of Womack et al. contains several cases where X(N2H+) is as low as 
or lower than 6.0 x 10""^^, e.g., p Oph and W49. Assuming a "canonical" CO abundance 




(8) 
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of 1 X 10~^, B95 also derived a value for X(N2H+) of 6.0 x 10"-^'^, more similar to the 
mean values of other cloud cores, by comparing column densities of N2H'^ with that of H2 
suggested by C^^O and C^^O data. B95 suggested this disparity could be reconciled if CO 
(and other molecules) in the non-outflowing gas of IRAS 4 were depleted relative to dust 
by factors of 10-20, but maintained abundances relative to each other that were similar to 
those seen elsewhere. 

Using Equation 8 and the above assumptions, the Mgas of 4A and 4B are 0.73 M0 and 
0.41 M0 respectively. As a comparison, recent SCUBA observations of continuum emission 
from dust around 4A and 4B by Sandell & Knee yield circumstcllar mass estimates of 
2.0-4.0 Mq and 0.65-1.3 M©, following their analysis but assuming a distance to NGC 1333 
of 350 pc and a range of Tk of 20-40 K. (In their analysis, Sandell & Knee assumed Tk = 
25 K and a 220 pc distance to derive masses.) The SCUBA-based masses are somewhat 
larger than ours, possibly due to their greater sensitivity to larger-scale emission from 
lower-density material. 

Another comparison of mass is provided by virial masses of the non-outflowing dense 
gas surrounding 4A and 4B. These can be calculated from the respective mean values of 
Ay and apparent sizes as 



for a uniform, self-gravitating sphere (e.g., sec Larson 1981). Using the mean values of 
listed in Table 4, and the values of R of 0.022 pc for 4A and 0.0093 pc for 4B (obtained 
from the extents of the integrated emission shown in Figure If) in Equation 9, virial 
masses of 1.1 Mq and 0.48 Mq are found for the non-outflowing dense gas around 4A and 
4B respectively, values quite similar to those determined above from the derived column 
densites. (Alternatively, values of X(N2H+) of 4.0 x 10~^^ and 5.1 x 10~^^, similar to that 
found by B95 for 4A, can be found for 4A and 4B respectively by inverting Equation 8 and 
assuming M^j^ = Mg^s-) 



Figures 7b and 7e show the distributions of Vlsr with position around 4A and 4B 
respectively obtained from hyperfine structure fitting. Systematic variations with position 
of the central velocity, VlsR) of the symmetric profiles of an optically-thin line can suggest 
velocity gradients in the gas. If a gradient is found symmetrically about an average velocity 




(9) 



4.3.2. Rotation 
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and position, it may be caused by rotational motions of the gas. 

As listed in Table 5, the mean Vlsr values of the non-outflowing dense gas associated 
with 4A and 4B are 6.95 km and 6.98 km s~^ respectively, each within the mean fitting 
uncertainty of the other (0.032 km s~"^). If these mean velocities of this dense gas are 
representative of the radial velocities of the 4A and 4B objects themselves, the 4A and 4B 
objects arc not moving relative to each other significantly along the line-of-sight, as might 
be expected if 4A and 4B were revolving about each other. (However, a Vlsr difference of 
0.28 km s^-*^ is found between 4A and 4B from modeling their respective associated inverse 
P-Cygni profiles; see §4.2.1, Table 3.) The average Vlsr obtained from the non-outflowing 
dense gas around both objects is 6.96 km s"^, which we take as the systemic velocity of 
IRAS 4. 

Near 4A, the rms variation of Vlsr is more than half the mean AV, i.e., 0.36 km s^^ vs. 
0.25 km s~^, suggesting the variations of Vlsr around 4A are signiflcant. Furthermore, the 
positional distribution of Vlsr in the non-outflowing dense gas about 4A shows systematic, 
symmetrical variations that may due to rotation. To flrst order, the southeastern half of 
this gas around 4A is blue-shifted (relative to the mean) by ~0.25-0.5 km s~^ and the 
northwestern half is red-shifted by ~0.0-0.5 km s~^, variations that could arise if this gas 
were rotating about 4A on an axis at a P.A. of ~45°. (Such an axis would be colinear with 
the "cavity" containing 4A and 4A'.) Rotating dense gas, however, cannot account for all 
variations of Vlsr about 4A. For example, two pockets to the north-by-northwest and one 
to the south-by-southwest have relative velocities in directions opposite to those expected 
from rotation at their locations. Given that the south-by-southwest pocket is unusually 
red-shifted and that it lies in the same direction of the blue 4 A outflow lobe (see §4.2.2), 
the pocket may be a local distortion of a rotational velocity held due to selective clearing 
of N2H+ at blue-shifted velocities by the outflow in that direction. A similar explanation 
may account for the unusually blue-shifted pockets in the same direction as the red outflow 
lobe. (The interpretation of velocity variations as caused by rotation is not unique, and the 
velocity distribution also could be due to bulk infall or outflow motions of gas on either side 
of 4A.) 

Near 4B, the characteristics of Vlsr are much different from those seen around 4A. The 
rms variation of Vlsr is less than half the mean AV^, i.e., 0.096 km s~^ vs. 0.25 km s^^, and 
only ~25% that of the non-outflowing dense gas around 4A, suggesting the variations of 
Vlsr are not significant. In addition, no systematic variation of Vlsr is seen with position, 
and hence no suggestion of rotation (and infall or outflow), at least from the N2H'*' 1-0 line. 
The difference in rotational properties between the non-outflowing dense gas surrounding 
4A and 4B is striking in the context of binary formation. Since 4A is a known binary of 1" 
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(or 350 AU) separation but 4B has not been confirmed to be a multiple, the relatively large 
rotation of the non-outfiowing dense gas about 4A suggests binaries may form preferentially 
in cores of relatively high intrinsic angular momentum but single stars form preferentially 
in cores of relatively low intrinsic angular momentum. (However, 4B may be yet a binary 
system, albeit with a much smaller separation than that of 4A and 4A'; see §4.2.2). 

4.3.3. Turbulent Motions 

NGC 1333 is located in the western Perseus molecular cloud, a region where a large 
number of stars are forming in close proximity from gas that is fairly turbulent, judging 
from the observed line widths. For example, IRAS 4A and 4B have among the largest 
N2H+ 1-0 line widths, 1.32 km s~^ and 1.04 km s^^ respectively, of the 47 Class or Class 
I objects within 400 pc sampled by M97. As the thermal FWHM of N2H+ is 0.247 km s"^ 
if TV = 40 K (i.e,. the maximum of the range derived by B95 for the non-outflowing 
gas), the broad lines seen by M97 in IRAS 4 suggest the presence of substantial turbulent 
motions in the dense gas on scales of 0.05 pc (i.e., the 27" IRAM 30-m Telescope beam 
FWHM at 93.17 GHz at 350 pc.) 

The PdBI data described here probes physical scales much smaller than those of 
single-dish observations, e.g., the synthesized beam FWHM is only 0.0055 x 0.0064 pc 
at 350 pc. Figures 7c and 7f show the distributions of the N2H+ 1-0 line width near 4A 
and 4B respectively, i.e., values of AV obtained from hyperfine structure fitting. Table 
5 lists the beam-averaged mean, rms, maximum and miniumum values of AV^ found for 
non-outflowing dense gas around 4A or 4B. From these data, a dramatic reduction in N2H"'" 
line width is seen, relative to the single-dish data. For example, maximum beam-averaged 
values of AV around both objects are <0.84 km s^^, and the mean beam-averaged values 
of AV are both only ~0.50 km s~^, much less than seen in single-dish data. 

At first glance, it is conceivable that the larger AV^ in single-dish observations could 
originate from the unresolved combination of intrinsically narrow lines with potentially 
wide variations of central line velocity in the dense gas (e.g., from rotation). However, the 
rms values of Vlsr and the mean values of AV are too small to support this possibility. 
Instead, the overall reduction of line width with scale in IRAS 4 may be a manifestation 
of the hnewidth-size relation, as observed in many molecular clouds (e.g., see McKee 1999 
and references therein), and the interferometer data is probing the smallest scales of this 
relation. At these scales, the nonthermal component of the line width is predicted to 
be comparable in magnitude to the thermal component, possibly due to dissipation of 
turbulence (Vazquez-Semadeni et al. 1999) or the cutoff of MHD wave propagation (Nakano 
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1998; Myers 1998) at small scales. 

A similar but less dramatic reduction of line width with scale was noted for 8 "quiescent 
cores" in the Serpens NW cluster region by WiUiams & Myers, which is presumed to be at 
a similar distance to NGC 1333 (i.e., 310 pc; de Lara, Chavarria-K, & Lopez-Molina 1991). 
In Serpens NW, the mean width ratio for lines observed between 10" and 50" resolution 

over all 8 cores was 0.59, whereas the mean width ratio for lines observed between ~4" and 
27" over 4A and 4B is 0.43, using the observations of M98 and the mean values of AV in 
Table 5. Interestingly, the 8 cores found by Williams & Myers have typical sizes of ~5000 
AU, line widths of ~0.5 km s^^, and infall speeds of '^0.5 km s^^, values all similar to those 
of the non-outfiowing dense gas around 4A and 4B (assuming the infall speeds found in 
§4.2.1 from H2CO apply to the N2H+ gas.) 

Despite the overall relative reduction in line width with scale seen here, the high spatial 
resolution of these observations reveals significant variations in AV" in the non-outfiowing 
dense gas around both objects (see Figures 7c and 7f), and a common narrow width of ~0.5 
km s"^ cannot describe the line emission. Figure 8 shows histograms of the distributions 
of AV as fractions of the total number of beam-averaged pixels around 4A (54 pixels or 
4B (12 pixels) in 0.15 km bins. Both have ranges of AV between 0.23 km s~^ to 0.83 
km s~^, with substantial fractions (15-30%) near 0.25 km s~^, the thermal line width at 
40 K, indicating the presence of very non-turbulent gas around each object. (Minimum 
beam-averaged AV values of 0.21 km s^^ near 4A and 0.16 km s~^ near 4B are found, 
located north of 4A and 4B, or at the peripherics of the dense gas.) Of the remainder, 
the distribution for 4A has peaks at 0.38 km s~^ and 0.68 km s~^, but that of 4B has 
only a peak at 0.23 km s~^. For 4A, hues of the first peak are located generally to the 
northwest and fines of the second peak are located generally to the southeast, although the 
populations are mixed. For 4B, the lines broaden systematically from north to south across 
its dense gas. 

The condensations from which 4A and 4B formed may be plausibly assumed to have 
had nearly-uniform initial velocity dispersions, as is seen in many starless dense cores (see 
Goodman et al. 1998; Caselh et al. 2001). The star formation process from that time until 
the present likely increased the velocity dispersion by differing amounts thoughout each 
condensation, due mainly to modest outfiow-core interaction or heating from the protostars, 
but very likely did not decrease the dispersion anywhere. Every independent measurement 
of line width in the gas around each object then gives an upper limit on the initial velocity 
dispersion of the original condensation. The typical upper limit comes from the mean 
present-day line width, ~0.5 km s~^ in both 4A and 4B, and the most constraining upper 
limit comes from the narrowest present-day line width, ~0.2 km s~^ in the gas around both 
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4 A and 4B. These new limits indicate that the initial level of turbulence in star-forming 

condensations may be very low, even in regions of clustered star formation. For example, if 
Tk = 20 K (the low end of from B95) and A\/(observed) = 0.225 km s~^ (the lowest 
AV bin in Figure 8), then (thermal) = 0.18 km s~^ for N2H"'" and Al^(nonthermal) = 
0.135 km s~^, as (nonthermal) = (AV^ (observed) - Al^^ (thermal)) 

If the velocity dispersion has evolved from its initial state though the injection of 
energy, line widths arguably should be largest at positions closest to possible sources of 
energy, and there is indeed some evidence for this in the dense gas around both objects. For 
example, there appears to be some spatial correspondence between large values of AV and 
the blue-shifted outflow features (as traced by H2CO) south of both 4A and 4B. However, 
the situation is not so clear cut and evidence to the contrary is also found. For example, the 
narrowest lines in the non-outflowing dense gas around 4A and 4B are found actually near 
the respective red outflow features north of each object. Moreover, large line widths are 
found across the southeastern half of the 4A dense gas, not just at positions near the blue 
outflow lobe. Given the factor of ~10-20 difference between the spectral width of the H2CO 
wings and the width of the N2H^ line even at positions of large AV^, whatever coupling 
does exist between the dense outflowing gas and N2H"'" appears to be very inefficient. 

The weak coupling between the fast outflow and the non-outflowing gas may be a 
consequence of their difference in speed, with greater coupling expected when the flow 
speed is more comparable to the velocity dispersion of its environment. On the other hand, 
possibly the N2H+ is depleted more rapidly than it can be accelerated to speeds similar to 
the outflow. (The contrast of widths cannot be ascribed to the outflow gas being less dense 
than the non-outflowing gas, as both transitions probe gas at n 10^~^ cm~^.) Similar 
contrasts in line width have been found toward other very young objects, but the proximity 
of outflow gas to the non-outflowing gas is not so clearly deflned in single-dish observations. 
For example, the survey of M97 contains several examples with outflow wings in H2CO and 
CS, but not in N2H+, on 0.02-0.05 pc scales, toward sources such as L483, S68N, or L1157 
(and, of course, 4A and 4B). The PdBl data here show emission from both outflowing and 
non-outflowing gas around 4A and 4B can be spatially coincident on scales of only ~0.006 
pc. 



5. Summary and Conclusions 

Millimeter interferometer observations are presented of NGC 1333 IRAS 4, a group 
of highly-embedded protostellar objects in the Perseus star-forming complex. Millimeter 
continuum data and line data allow the immediate dense gas and dust environments of 
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its protostellar objects to be characterized simultaneously at high spatial resolution. The 
major results found from these observations include: 

1. Inverse P-Cygni profiles of H2CO 3i2-2ii that show red-shifted absorption and 
blue-shifted emission have been detected toward 4 A and 4B. Inverse P-Cygni profiles of 
CS 3-2 and N2H+ 1-0 also have been detected toward 4A. These profiles provide the least 
ambiguous evidence for infalling gas motions associated with low-mass young stellar objects. 
The depth of the absorption in the H2CO line are equal to the A = 1.3 mm continuum 
intensities of the protostars along the same lincs-of-sight. Modeling the H2CO profiles with 
a "two-layer" radiative transfer code, infall velocities toward 4A and 4B of 0.68 km 

and 0.47 km s^^ are estimated respectively. From these, masses of 0.53 Mq and 0.17 Mq, 
and accretion rates, Mi^, of 1.1 x 10~^ Mq yr~^ and 3.7 x 10~^ Mq yr~^ for 4A and 4B 
are estimated respectively, given assumptions abut ri„ and n. Accretion timescales for both 
objects are similar, with tacc — M/M of 4700 yr and 4600 yr for 4A and 4B respectively. 

2. Outfiowing gas from both 4A and 4B is detected from bright wing emission of H2CO 
3i2-2ii and CS 3-2. The outfiows from each source differ, in that only the southern blue 
lobe from the extended outfiow from 4 A is detected but 4-5 compact features from 4B 
are detected. For 4A, the velocity structure of the emission suggests the H2CO is tracing 
outfiow cavity walls. The orientation of the cavity walls suggests the outfiow has precessed 
from orientations evident at larger scales, due possibly to the binarity of 4A. For 4B, the 
multipolar structure could be also from cavity walls of an outfiow that is more compact 
than that from 4A. However, the complex structure may originate from several juxtaposed 
lobes from a multiple object. 

3. Non-outfiowing dense gas surrounding 4A and 4B individually is detected from 
N2H+ 1-0 emission. With assumptions, we determine mean column densities of 2.4 x 10^^ 
cm~^ and 3.9 x 10"*^^ cm~^ respectively. Furthermore, masses of this gas arc estimated to 
be 0.73 Mq and 0.41 Mq respectively, comparable to other estimates of circumstcUar mass 
derived from submillimeter dust emission and the virial theorem. The variations of Vlsr 
around each object suggest the non-outfiowing dense gas about 4A is rotating, but the data 
do not suggest similar motions about 4B. 

4. The widths of the N2H+ 1-0 fine in the non-outfiowing dense gas around both 4A 
and 4B are smaller by factors of 2-3 on average than those observed in single-dish data, 
due to the reduction of turbulent motions on smaller physical size scales. The data reveal 
a new, low upper bound to the initial velocity dispersion in the protostellar core from 
which IRAS 4 formed. For example, some regions around each object (~15-30% of the 
projected surface areas of line emission) exhibit lines with very narrow widths, e.g., ~0.2 
km s~^. Such narrow line widths are consistent with the incoherent motions of the gas 
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having a nonthermal component of a magnitude comparable to that of thermal broadening 

at expected gas temperatures, and may represent the initial state of the core gas prior to 
the formation of the protostars or their outflows. The lines show considerable variation 
of width, however, even on the small scales probed. Coupling between outflow and the 
non-outflowing dense gas, that could inject turbulent energy into any initially non-turbulent 
gas, appears highly inefficient, as line widths from H2CO and N2H+ can differ by factors of 
10-20 on scales of only 0.006 pc. 
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A. High-Resolution Continuum Observations of NGC 1333 IRAS 4 

Table 2 lists the flux densities, peak intensities, and deconvolved sizes of each source, 
derived from fltting single-component two-dimensional Gaussians to observations of each 
source individually at each wavelength with the AlPS task IMFIT. (Positions are also 
deflned by these flts, e.g., those shown in Figure 1 are those of the peak continuum intensties 
at A = 1.3 mm from Gaussian fitting.) Flux density uncertainties are comprised of those 
derived from the fits, an assumed 15% uncertainty in the absolute calibration at each 
wavelength, and an estimated 10% measurement uncertainty. Values in Table 2 for 4B' at A 
= 1.3 mm are particularly uncertain due to its location at the mosaic edge (see Figure la), 
although the peak position obtained from the fit agrees quite well with those from other 
wavelengths. 

Figure Al shows the millimeter wavelength spectral energy distributions (SEDs) of 
each IRAS 4 object from our continuum data. In addition. Figure Al shows data obtained 
at A = 800 fim by Lay et al. from the JCMT-CSO Interferometer, at A = 1.3 cm and A 
= 3.6 cm by Mundy et al. (1993) at the VLA, and at A = 3.6 cm and A = 6.0 cm by 
Rodriguez et al. (1999) at the VLA. At all wavelengths, 4A is brighter than 4B, which is in 
turn brighter than 4B'. The SEDs of Figure Al are reasonably well fitted by straight lines 
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of similar slopes, i.e.,. = -2.7 ± 0.05 for 4A, -2.9 ± 0.14 for 4B, and -2.6 ± 0.23 for 4B', 
using linear least-squares fitting over the range 1.3 mm < A < 1.3 cm. (For 4A and 4B, 
only our data at A = 1.3 cm was included for the fit, as they are similar (i.e., within errors) 
to those found by Mundy et al. at the same wavelength.) For 4B', the linear fit included 
the fiux density at A = 1.3 mm, although this measurement is particularly uncertain, and 
excluded our upper limit at A = 1.3 cm. 

The similar slopes reveal that continuum emission for the IRAS 4 objects at A < 1.3 
cm is dominated by thermal emission from dust grains surrounding each source. Given 
the deconvolved sizes of the emission, the warm dust being traced is likely situated within 
accretion disks or the inner radii of envelopes around each source. If the dust emission is 
optically thin (and we assume the same region is probed) at all wavelengths, the slope, a, 
can be used to estimate the indices of the mass opacity of dust, i.e., (3. As -a — (3 -\- 2 
in the Ray leigh- Jeans limit, dust opacity indices of 0.7 ± 0.05 for 4A, 0.9 ± 0.14 for 4B, 
and 0.6 ± 0.25 for 4B' are found from our ai, each consistent with a common value of (3 
~ 0.75, and significantly less than the (3 = 2 value assumed typically for dust in the ISM. 
These values are similar to the (3^1 values inferred for dust around many other young 
stellar objects, a decrease in (3 possibly due to growth or composition changes of the dust 
(Beckwith & Sargent 1991). 

The submillimeter wavelength data of Lay et al., not included in the linear fits shown 
in Figure Al, deviate from the derived fits because it was obtained at a resolution of 
~1", and much of the emission from dust in the inner radii of envelopes may have been 
resolved out in those data. The centimeter wavelength data of Rodriguez et al., also not 
included in the linear fits, deviate dramatically from the fines fit solely to the data at A < 
1.3 cm, suggesting that free-free emission becomes increasingly dominant in the SEDs of 
these sources at long wavelengths. Extrapolating these data to shorter wavelengths, the 
contribution of this emission to the respective flux densities is small, i.e., <20% at A = 1.3 
cm and <5% at A = 6.9 mm for both 4A or 4B. 

Departures of individual flux densities at millimeter wavelengths from the linear flts 
are seen least in the flt to the 4B SED, whose reduced — 1-6, but are seen more for 
the flts to the 4A SED (reduced = 7.9) and 4B' SED (reduced = 9.5). These latter 
large values of reduced be due to underestimates of the absolute flux calibration 

uncertainties. 

At speciflc wavelengths, the largest departures from the linear flts are the positive ones 
found at A = 3.2 mm for all 3 objects. The object 3C 84 was also observed during all tracks 
as IRAS 4, but a similar positive departure of its fiux density at A = 3.2 from a linear fit 
is not seen. This consistency suggests the deviations in the IRAS 4 SEDs from their fits at 
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this wavelength may be real, and the slopes of the IRAS 4 SEDs at A < 3.2 mm may be 
shallower than those found in earlier fits. Dust emission at A > 3.2 mm may be declining 
at a steeper overall slope, possibly because the cooler, outer layers of the IRAS 4 envelopes 
are resolved out at our high resolutions, or have lower optical depths than the inner regions. 

Slightly shallower slopes of the linear fits are found for the 4A and 4B SEDs if our 
VLA data are additionally excluded from the fits, i.e., = -2.6 ± 0.29 and -2.7 ± 0.30 
respectively, but these do not differ significantly from slopes obtained with the VLA data 
included. One slope that does differ significantly is that of 4B', namely a2 = -1.8 ± 0.39. 
The slopes of all three SEDs change only slightly if the NMA data arc also excluded. 
Unlike for 4A and 4B, both ai and a2 for 4B' are significantly more than the values of 
-3.8 or -3.9 derived by Choi (2001) from interferometer continuum data at A = 3.4 mm 
and A = 2.7 mm. Values of /3 derived from our slopes suggest so 4B' could have dust with 
significantly different opacity than 4 A and 4B, possibly from a population of larger dust 
grains. Alternatively, the dust may be optically thick. Given the relative dimness of 4B' 
relative to 4A or 4B, more sensitive data than those shown here (especially at A = 1.3 mm 
and A = 6.9 mm) are needed to properly evaluate either possibility. 
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Table 1: Log of Observations 



A V Interferometer Configuration Observing Dates Beam FWHM Continuum rms 

^mm) (GHz) (") (in.Ty l)eain--^) 



1.3 225.7 IRAM PdBI 



5C2-W09 
4ant- special 

5D-N09 
4D2+N13 



21 Jul 1997 
19 Aug 1997 
19 Sep 1997 
14 Nov 1997 



2.0x1.7 



17.5 



2.0 147.0 



NMA 



D 

C 



27 Dec 1997 
09 Mar 1998 



3.2x2.3^ 



14.0 



2.2 137.0 



NMA 



D 

C 



27 Dec 1997 
09 Mar 1998 



3.4x2.6 



10.5 



3.2 



93.2 



IRAM PdBI 



5C2-W09 
5C2-W09 
4ant-special 

5D-N09 
4D2+N13 



20 Jul 1997 

21 Jul 1997 
19 Aug 1997 
19 Sep 1997 
14 Nov 1997 



2.9x2.6 



4.5 



6.9 



43.3 



13.3 22.5 



VLA 



VLA 



D 



D 



25 Nov 1997 



25 Nov 1997 



2.1x1.9 
3.8x3.1 



0.30 
0.28 



"CS 3-2 data tapered to a 6'.'5 x 4'.'4 FWHM beam for Figure le and analysis. 
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Table 2: Integrated Flux Densities, Peak Intensities and Deconvolved Sizes 











Deconvolved Size 


Object 


A 


i^total 


Tpeak 

V 


Major X Minor Axis 


P. A. 




(mm) 


(mJy) 


(mJy beam~^) 


/ A T T At t\ 

(AU X AU) 


(degrees) 


4A 


1.3 


3100 ± 470 


1200 ± 16. 


920 X 720 


132 




2.0 


690 ± 110 


450 ± 14. 


770 X 590 


129 




2.2 


580 ± 90. 


420 ± 10. 


740 X 570 


119 




3.2 


320 ± 49. 


200 ± 4.4 


790 X 650 


177 




6.9 


26. ± 4.0 


16. ± 0.29 


580 X 490 


138 




13.3 


4.7 ± 0.95 


3.1 ± 0.27 


1100 X 630 


125 


4B 


1.3 


1100 ± 170 


590 ± 17. 


600 X 560 


151 




2.0 


320 ± 56. 


250 ± 14. 


630 X 390 


119 




2.2 


240 ± 41. 


210 ± 11. 


440 X 300 


151 




3.2 


93. ± 16. 


90. ± 4.5 


<570 X <570 






6.9 


7.8 ± 1.3 


5.8 ± 0.29 


630 X 110 


128 




13.3 


1.2 ± 0.56 


1.0 ± 0.28 


<1400 X <1400 




4B' 


1.3 


230 ± 47. 


200 ± 18. 


<640 X <640 






2.0 


79. ± 29. 


70. ± 14. 


<840 X <840 






2.2 


82. ± 21. 


89. ± 10. 


<780 X <780 






3.2 


45. ± 11. 


41. ± 4.5 


<690 X <690 






6.9 


2.2 ± 0.65 


1.9 ± 0.30 


<840 X <840 






13.3 


<0.84 


<0.84 
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Table 3: "Two-layer" Model Parameters for Inverse P-Cygni Profiles^ 



Object To ^ Jc Jf Jr VlSR Cr Vin 

(K) (K) (K) (kms-^) (km s"^) (km s"^) 



4A 2.0 0.3 45 3 19 6.85 0.34 0.68 

4B 1.6 0.3 30 3 20 7.13 0.51 0.47 



"see Figure 4 for model profiles. 
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Table 4: Physical Quantities of Mass Flows in NGC 1333 IRAS 4 



Derived 






IRAS 


IRAS 


Quantity 


Abbreviation 


(units) 


4A^ 


4B^ 


InfalP: 


Mass Accretion Rate 




M© yr-1 


1.1 X 10"^ 


3.7 X 10-5 


Outflou^: 


Mass 


Mout 


Mo 


6.9 X 10-^ 


1.5 X 10"^ 


Momentum 


p 

out 


Mq km s~^ 


1.6 X 10-2 


1.6 X 10-^ 


Kinetic Energy 


Eout 


erg 


4.9 X 10^2 


3.2 X 10^1 


Momentum Flux 


FH2C0 


Mq km s^^ jY^^ 


4.2 X 10-5 


3.2 X 10-<^ 


Mechanical Luminosity 


-^mech 




1.3 X 10-1 


7.4 X 10-^ 



"southern blue lobe only. 

''includes contribution from eastern blue lobe, some of which may originate from 4A'; see text, 
'^derived from very uncertain ri„; see text. 

''for dense gas in outflow only; assumes outflow inclination angle of 10° for both 4A and 4B. 
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Table 5: N2H+ Line Properties Derived from Hyperfine Structure Fits 



Value 


Unit 


mean 


4A 
rms 


min 


max 


mean 


4B 
rms 


min 


max 




K 


7.8 


2.1 


3.6 


13. 


12. 


4.4 


4.8 


21. 


Vlsr 


km 


6.95 


0.359 


6.48 


8.2 


6.98 


0.0928 


6.75 


7.38 




km 
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Fig. 1. — Millimeter continuum and line interferometer maps of NGC 1333 IRAS 4. Positions 

of contimiTim somxcs 4A, 4B, and 4B' at A = 1.3 mm are denoted by large black stars while 
the position of continuum source 4A' from Looney et al. (2000) is denoted by a smaller 
black star. The black ellipses placed in the bottom left of each panel represent the FWHMs 
of the respective synthesized beams. Outer dashed contours represent the masking cutoffs 
at a specific gain level of the mosaics: 20% for the IRAM PdBI, or 50% for the NRO NMA 
data, a) X — 1.3 mm continuum emission. Plain contours start at 52.5 mJy beam~^ (3 
a), increasing in steps of 2 cr to 297.5 mJy beam~^; the dashed contour denotes -52.5 mJy 
beam~^. b) X — 2.0 mm continuum emission. Plain contours start at 42.0 mJy beam~^ (3 
a), increasing in steps of 2 cr to 126.0 mJy beam~^; the dashed contour denotes -42.0 mJy 
beam~^. c) X — 3.2 mm continuum emission. Plain contours start at 13.5 mJy beam~^ 
(3 a), increasing in steps of 2 o" to 49.5 mJy beam^^; dashed contour denotes -13.5 mJy 
beam^^. The arrows to the northeast and southwest of 4A+4A' denote the directions of the 
large-scale CO outflow detected by Blake et al. (1995) d) The integrated intensity of H2CO 
3i2^2ii emission ( "zeroth-moment" ; see text). The bold contour shows the half-maximum 
isophote (57 K km s~^), while plain contours step up and down in increments of 1.5 a (10.5 
K km s~^) from this level, e) The integrated intensity of CS 3-2 emission. The bold contour 
shows the half-maximum isophote (11.1 K km s~^), while plain contours step up and down 
in increments of 1.5 a (4.5 K km s~^) from this level. /) The integrated intensity of N2H"'" 
1-0 emission (over all 7 hyperfine components). The bold contour shows the half-maximum 
isophote (7.6 K km s^^), while plain contours step up and down in increments of 1.5 a (1.2 
K km s~^) from this level. 

Fig. 2. — Grids of H2CO 3i2-2ii spectra at selected positions at and around 4A and 4B. In 
both grids, panels are spaced by 2" in R.A. or Dec, ~1 FWHM of the beam at 225.7 GHz, 
the H2CO 3i2-2ii rest frequency, so adjacent panels in the H2CO grid are independent. The 
dashed line represents the determined Vlsr of the IRAS 4 system of 6.96 km s~^. a) top 
grid: Panels centered at the 4A position at A = 1.3 mm (R.A. (2000) = 03:29:10.5, Dec. 
(2000) = +31:13:31.4). b) bottom grid: Panels centered at the 4B position of 4B at A = 1.3 
mm (R.A. (2000) = 03:29:12.0, Dec. (2000) = +31:13:08.2). 

Fig. 3. — Grids of N2H+ 1-0 spectra at the same positions at and around 4A and 4B defined 
for Figure 2. Due to the larger beam at 93.17 GHz, the N2H+ 1-0 rest frequency, the adjacent 
panels are not independent, as in Figure 2. In addition, the velocity range in each panel 
here is wider relative that of to Figure 2, although panels in both Figures are centered at 
the same velocity. 

Fig. 4. — H2CO 3i2-2ii spectra from the IRAM PdBI toward the positions of peak continuum 
intensity near 4A [panel a) and 4B [panel b). The dotted line denotes the average systemic 
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Vlsr of +6.96 km s ^. In each panel, a dashed line shows a spectrum resulting from a 
simple two-layer model of infall toward each object. Respective model parameters are listed 
in Table 3. 

Fig. 5. — Channel maps of H2CO 3i2-2ii emission revealing the underlying structure of the 
southern blue lobe of the 4A dense gas outflow. Each channel is the mean of 3 channels of the 
original data cube, and has a channel width of 0.3113 km s~^. The velocity of a given panel 
is shown at the upper left in km s~^; only a fraction of channels blueward of the systemic 
velocity are shown. Continuum positions of 4A and 4A' are denoted by large stars and small 
stars respectively. Solid contours begin at 3.5 K and increase in steps of 3.5 K, and dashed 
contours begin at -3.5 K and decrease in steps of -3.5 K. The grey scale ranges from -3.5 K 
to 28.2 K. 

Fig. 6. — Channel maps of H2CO 3i2-2ii emission depicting the red and blue lobes of the 
4B dense gas outflow. Channel width, contours, grey scales, and symbols are the same as 
described for Figure 5, except stars denote the continuum locations of 4B (near center) and 
4B' (to the left). In addition, an equal number of channels redward and blueward of the 
systemic velocity of 4B are shown. 

Fig. 7.— The spatial distribution of T|f*, Vlsr and AV of the N2H+ 1-0 line around IRAS 4A 
and 4B, obtained by fitting the 7 hyperfine components of the line simultaneously with the 
HFS fitting routine of CLASS. Only positions where the peak intensity of brightest hyperfine 
component was >5x the spectral baseline rms are shown, and data at other positions have 
been blanked. The positions of continuum objects 4A, 4A', and 4B, are represented by stars, 
and are the same positions as denoted in Figure 1. a) Distribution of Tg"* around 4A, where 
the color scale is defined from K to 40 K, and contours begin at 2.5 K and continue in 
increments of 2.5 K to 40.0 K. b) Distribution of Vlsr around 4A, where the color scale is 
defined from 6.0 km s~^ to 8.0 km s~^, and contours begin at 6.2 km s~^ and increase in 
increments of 0.2 km s~^ to 7.8 km s~^. c) Distribution of AV around 4A, where the color 
scale is defined from 0.1 km s^^ to 1.5 km s^^, and contours begin at 0.2 km s~^ and increase 
in increments of 0.2 km s~^ to 1.4 km s~^. d) Distribution of T^°* around 4B, with the same 
color scale and contour levels defined for panel a), e) Distribution of Vlsr around 4B, with 
the same color scale and contour levels defined for panel b). /) Distribution of AV around 
4B, with the same color scale and contour levels defined for panel c). 

Fig. 8. — Distribution of AV from N2H+ 1-0 emission around 4A (panel a) and 4B [panel 
b) in 0.15 km s~^ wide bins as fractions of the total number of non-blanked pixels, i.e., 54 
pixels for 4A and 12 pixels for 4B. (The criterion for blanking is the same as for Figure 7.) 
Values of AV are obtained by deconvolving the measured values by the 0.16 km s~^ velocity 
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resolution of the PdBI correlator. The dashed line indicates the FWHM expected from the 
N2H+ 1-0 line from thermal broadening only at = 40 K, the maximum estimated 
for the non-outflowing gas around 4A by B95. The dot-dashed line indicates the mean AV 
measured from spectra around each protostellar object. 

Fig. Al. — Spectral energy distributions of 4A, 4B, and 4B' objects from interferometer 

continuum observations between A = 800 /xm to A = 6 cm. Filled squares, filled triangles, 
or filled stars show data obtained respectively from the IRAM PdBI, the NRO NMA, or 
the NRAO VLA for this paper. The open circles show data at A = 800 fim from the 
JCMT-CSO Interferometer by Lay, Carlstrom, & Hills (1995) for 4A and 4B. Open squares 
show data at A = 1.3 cm and A = 3.6 cm from Mundy et al. (1993) and open diamonds 
show data at A = 3.6 cm and A = 6.0 cm by Rodriguez et al. (1999) for 4A and 4B all from 
the Very Large Array. Downward-pointing arrows indicate 3 a upper limits. Error bars are 
1 (7 (see text and Table 2). Solid lines indicate lines of slope ai from resulting from linear 
least-squares fits to our flux densities of each object at 1.3 mm < A < 1.3 cm. Dashed lines 
indicate lines of slope a2 resulting from similar fitting but using only our data at 1.3 mm < 
A < 3.2 mm. 
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